Out-of-plane nesting driven spin spiral in ultrathin Fe/Cu(001) films 
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Epitaxial ultrathin Fe films on fee Cu(OOl) exhibit a spin spiral (SS), in contrast to the ferromagnetism of 
bulk bee Fe. We study the in-plane and out-of-plane Fermi surfaces (FSs) of the SS in 8 monolayer Fe/Cu(001) 
films using energy dependent soft x-ray momentum-resolved photoemission spectroscopy. We show that the SS 
originates in nested regions confined to out-of-plane FSs, which are drastically modified compared to in-plane 
FSs. From precise reciprocal space maps in successive zones, we obtain the associated real space compressive 
strain of 1.5±0.5% along c-axis. An autocorrelation analysis quantifies the incommensurate ordering vector 
g=(27r/a)(0,0,^0.86), favoring a SS and consistent with magneto-optic Kerr effect experiments. The results 
reveal the importance of in-plane and out-of-plane FS mapping for ultrathin films. 



Fe metal is delicately perched on a magneto- structural in- 
stability, with antiferromagnetic Cr and Mn to its left, and fer- 
romagnetic Co and Ni to its right in the periodic table. In spite 
of its high ferromagnetic Curie temperature (Xc=1043 K) and 
large magnetic moment (=2.22 //b) in its bulk body-centered- 
cubic (bec) form, theoretical studies have predicted face- 
centered-cubic (fee) Fe to be nonmagnetic, ferromagnetic, an- 
tiferromagnetic or exhibit a spin spiral (SS: corresponding to 
a periodic variation of the angle of the spins), sensitively de- 
pending on its lattice parameter (TJ [21 |4). The phase di- 
agram of Fe metal is indeed very rich with, (1) the bec to 
fee Bain transformation at 1184 K (5), (2) a hexagonal-close- 
packed (hep) non-magnetic phase at high pressure (P) and low 
temperature (T) which shows superconductivity [6], while (3) 
the high-P, high-T phase was debated in terms of a hep or 
orthorhombic structure [7]. Also, fee Fe nanoparticle precip- 
itates (diameter ~ 50 nm) stabilized in a Cu matrix shows a 
SS state with an ordering vector q = (27r/a)(0. 12,0, 1.0) and 
a low-spin moment (=0.7 /ie) per atom [ 8] [9|. Remarkable 
progress has been achieved in understanding the interplay of 
crystal structure and electronic instabilities in realizing spin 
order and its influence on properties of complex materials 
such as multiferroic oxides, superconducting copper-oxides 
ifTOl and iron-pnictides [11], etc. Surprisingly, despite exten- 
sive efforts, the simpler SS associated with the bcc-fcc (a-7) 
transition in elemental Fe and its momentum-resolved elec- 
tronic structure has remained a challenging unsolved problem 
for over 40 years (TJ. 

In contrast to experiments at high-T, high-P or nanopar- 
ticles confined in a matrix, a simpler route was developed 
to study structure property correlations in Fe: Fe thin films 
grown epitaxially on fee Cu(001) substrate fl2l [T3l fT4l [Bl . 
Fe thin films on Cu(001) also show a complex magnetic and 
structural phase diagram (T21 031 021 021 03) and it is known 
that (1) Below 4 monolayers (MLs), it is ferromagnetic and 
has the face-centered- tetragonal (fct) structure (Region I). (2) 
Between 5-11 ML, it has the fee structure and a SS, with a top 
bilayer which is ferromagnetic (Region II). The SS ordering T 
is Tss-200 K CUCS). ( 3 ) Above 12 ML > il transforms to the 
bulk ferromagnetic bec structure (Region III). The SS phase 



of Region II is the focus in this work, and as determined by 
magneto-optic Kerr effect measurements, it has a SS ordering 
vector of q z = 2nl2.6d or 2i\l27ld (where d is the interlayer 
distance of the Fe thin film) ||T6j |T8]|. This corresponds to q 
= (27r/a)(0,0,^0.75) (where a = 2d is the lattice parameter of 
the Fe thin films). The magneto-optic Kerr effect (MOKE) ex- 
periments provided the z(out-of-plane)-component of the SS 
vector and ruled out a type-1 antiferromagnetic structure. The 
obtained average magnetic moment per atom was ~1.5 //b> 
significantly larger than for the nanoprecipitates [021 [18]]. A 
depth resolved X-ray magnetic circular dichroism study indi- 
cated a magnetic moment of 1 .7 [1b for the SS phase, thus con- 
firming a clear difference compared to the nanoprecipitates 
fl9l . The possibility of a SS in fee Fe has been extensively 
discussed from theory. The salient features for fee Fe from 
theory indicates absence of Fermi surface (FS) nesting but al- 
lows for SSs with ordering vectors qi = (27r/a)(0,0,0.6) or q 2 
= (27r/a)(0.5,0,1.0) EH. Both these SSs show high- spin mo- 
ment (~1.5 /ie) and total energy considerations indicate qi is 
slightly more stable than q 2 . Also, a spin density wave (SDW: 
corresponding to a periodic modulation in the magnitude of 
collinear spins) was found to be extremely unfavourable com- 
pared to a non-collinear SS [20]. It is noted that qi and qr 2 
do not match with q of the MOKE measurements. The most 
serious limitation to date, however, is the absence of experi- 
mental results relating the momentum (k) resolved electronic 
structure to the SS in epitaxial ultrathin Fe/Cu(001) films. 

In order to make this connection, it is necessary to 
determine the in-plane and out-of-plane FSs of ultrathin 
Fe/Cu(001) films. While in-plane FSs are routinely measured 
with a fixed photon energy and angle-resolved photoemis- 
sion spectroscopy (ARPES), the specific requirement of in- 
vestigating the q z component of the SS in the out-of-plane 
FSs necessarily requires ARPES with tunable energy I2TH221 . 
While crystalline order along the c-axis has been reported for 
6-8 ML Fe/Cu(001) films using surface diffraction Q21, the 
present experiments constitute the first out-of-plane FS mea- 
surements for ultrathin films. Early study of the thickness de- 
pendent valence band structure of Fe/Cu(001) films identified 
distinct changes at 5 ML thickness, the borderline of regions 
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I and II l23ll . FS mapping carried out with a fixed energy (hv 
= 75 eV), showed diffuse in-plane FSs which could not be 
satisfactorily explained by generalized gradient approxima- 
tion band structure calculations of the bulk bcc or fee struc- 
ture. The added complication of the top ferromagnetic bi- 
layer motivated our use of soft x-ray (SX)-ARPES as the tech- 
nique of choice due to its larger probing depth, typically 10- 
20 A. Recent studies using SX-ARPES have reported the 3- 
dimensional (D) momentum resolved electronic structure of 
solids EH [22l [241 • In a very recent work, it was used to ob- 
tain the 3D charge density of bulk copper metal l25l . Hence, 
in comparison to low-energy ARPES, energy-dependent SX- 
ARPES (hv ~ 400-800 eV) is well-suited to investigate the 
3D FSs of the SS phase in ultrathin Fe/Cu(001) films. 
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FIG. 1 . Characterization of the sample and determination of k z . (a) 
Typical RHEED oscillations for Fe/Cu(001) grown at RT. (b) and (c) 
LEED patterns for clean Cu(001) substrate and 8 ML Fe thin films 
on Cu(001). (d) The volume BZ of fee Fe. The in-plane and out-of- 
plane regions probed in the present study are marked by blue and red 
lines, respectively, (e) EDCs along in-plane T-X (k x ) measured at 
hv = 430 eV. (f) and (g) EDCs along out-of-plane F-X (k z ) measured 
using hv = 415-580 and 740-940 eV. (h) The calculated k z positions 
from the photon energy dependence of EDCs. We used an inner 
potential and work function of 6 and 5 eV, respectively. The colors 
of the arcs corresponds to EDCs in (f) and (g). The EDCs plotted as 
thick lines correspond to momenta associated with the SS, as marked 
by arrows in Fig.[3|e). 



Epitaxial Fe films were grown on Cu(001) using electron 
beam deposition at room temperature (RT) (26). The thick- 
ness of the Fe thin films was controlled to 8 ML by mon- 
itoring reflection high energy electron diffraction (RHEED) 
oscillations (Fig. [IJa)). The crystallinity and crystal orienta- 
tion of the Cu(001) substrate and the deposited Fe thin films 
were measured by low energy electron diffraction (LEED) 



(Fig. [TJb) and (c), respectively). The RHEED oscillations 
and LEED patterns are in very good agreement with earlier 
work Q31Q21IT8]], indicating the high quality of the epitaxial 
ultrathin 8 ML Fe/Cu(001) films. 
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FIG. 2. ARPES results of Fe(8 ML)/Cu(001) at 50 K. (a) and (d) In- 
plane band maps measured at hv = 430 and 580 eV. The gray lines 
and points are the band calculation of fee Fe for q = (27r/a)(0.5,0,l) 
|4|. (b) and (e) EDCs at some k x points for each photon energy. The 
numbers labeled on EDCs indicate the probed k x positions, as shown 
in (a) and (d). (c) and (f) MDCs at E F . 



SX-ARPES measurements were performed with an elec- 
tron energy analyzer, Gammadata-Scienta SES2002, at undu- 
lator beamline BL17SU of SPring-8 using a grazing incidence 
geometry (<5°) spectrometer (26j [27 , 28]. The grazing inci- 
dence spectrometer makes the ARPES measurements highly 
efficient and also ensures that the x-ray photon momentum 
(~0.1T-X at hv = 400 eV) imparted to the electrons are ac- 
counted for easily |22l[24]|. Circularly polarized x-rays were 
used to avoid the symmetry selectivity of linearly polarized 
x-rays l26l . The measurements were carried out at 50 K, 
which is well below Tss~200 K, to minimize indirect transi- 
tion losses. The k x and k z positions in the 3D Brillouin zone 
(BZ) (Fig.[TJd)) were accurately determined from energy and 
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momentum distribution curves (EDCs and MDCs) obtained 
from intensity maps measured for incident photon energies 
probing two successive BZs. The arcs in k x -k z momentum 
space traced by the specific photon energies used are shown 
in Fig.[TJh). Figure[TJe) shows the EDCs for the k x scan along 
T-X at hv = 430 eV, while Fig. [iff) and (g) show EDCs for 
T-X along k z in successive BZs, respectively, as marked in 
Fig.[TJh). The overall similarity of the EDCs, particularly at 
the high symmetry points T and X along k x and k z confirms 
the viability of 3D FS mapping of the Fe films. However, since 
we are measuring ultrathin films with soft x-rays, we have to 
be careful about contributions from the copper substrate in the 
EDCs. It is well-known that copper shows very weak inten- 
sity due to a s band crossing near X point (along T-X and 
X-W) within a binding energy of 2 eV [26, 29 ]. We have ac- 
tually confirmed the same by measuring the clean copper sub- 
strate ARPES spectra before the film depositions l26l . Then, 
as we show in the following results, the SS associated elec- 
tronic states are well separated from the copper states in en- 
ergy and /c-space. Thus, the intensity maps within 2 eV bind- 
ing energy are dominated by the Fe spectral weight. Fig.|2ja) 
and (b) show the in-plane T-X and X-W intensity maps of the 
Fe films, obtained using hv = 430 eV and 580 eV, respec- 
tively. Figure |2jb)/(e) and (c)/(f) are the corresponding EDCs 
and MDCs which show clear band crossings. Early work as- 
signed the band crossings along T-X to the Ai and A 5 bands, 
as suggested by calculations for the parent high- spin fee Fe 
(30). A comparison with calculated band dispersions for the 
SS vectors obtained in the local spin density approximation 
(LSDA), however, does not give a satisfactory match with the 
data for q 1 = (27r/a)(0,0,0.6). For q 2 = (27r/a)(0.5,0,1.0), the 
calculated data are overlaid in Fig. [2ja) and (b) after shift- 
ing the energy by 0.4 eV toward Fermi edge. While the data 
show similarities with calculations, e.g., MDCs show band 
crossings matching the calculation for the X-W cut or for the 
crossing at 0.35 A -1 along T-X, there are clear discrepancies 
also observed: the intense band crossing at 1.05 A -1 along 
T-X and the band dispersions between 0.5 to 2 eV binding en- 
ergies do not match well with calculations. Thus, while band 
dispersions and crossings can be clearly observed, the experi- 
mental results do not match the LSDA results for the SSs with 
q x = (27r/a)(0,0,0.6) and q 2 = (27r/a)(0.5,0,1.0). 

We then embarked on a determination of the in-plane (k x - 
k y ) and out-of-plane (k x -k z ) FSs of the Fe films in order to 
trace the electronic structure associated with the SS. The in- 
plane (Fig.|3ja)) and out-of-plane (Fig.[3fb)) intensity maps at 
E F were obtained from polar angle scans, at hv = 430 eV and 
as a function of energy hv = 385-595 eV, respectively. The 
raw data was measured over more than one-quarter of the BZ 
(Fig.[3ja) and (b)) [26]. From peaks in MDCs, we have deter- 
mined the FS crossings and symmetrized to obtain the exper- 
imental in-plane and out-of-plane FSs for full BZ (Fig. [3jc) 
and (d)). Since the fee structure is equivalent along (100), 
(010) and (001), the in-plane and out-of-plane FSs should be 
equivalent in the absence of electronic/magnetic instabilities 
along preferred directions. The results instead show remark- 



able differences in the in-plane and out-of-plane FSs. The in- 
plane FSs show a four- fold symmetry while the out-of-plane 
FSs show two-fold symmetry. While the in-plane FSs show no 
evidence for nesting, a clear nesting is observed in the out-of- 
plane FSs. This is directly seen in the MDC plot comparison 
for the k x , k y and k z cuts along T-X (Fig. J3je)). While the 
k x and k y cuts are very similar, the k z cut shows one addi- 
tional peak between T and X. Thus, the peak-to-peak separa- 
tion for the symmetrized data, marked by arrows in Fig.|3je), 
corresponds to (27r/a)x(0.86±0.1). The additional peaks in 
the MDCs are observed over an extended region in k x (red 
dashed lines in Fig. [3jd)), thus identifying the q z component 
of the nesting vector. The identification of the SS nesting 
vector in the out-of-plane FSs is based on a precise determi- 
nation of reciprocal-space maps in successive BZs. In turn, 
this provides evidence for a subtle structural modification in 
real space, namely, a small but finite compressive strain of 
1.5±0.5 % along the c-axis l26l . 

In order to determine the q x component of the SS vector, we 
carried out an autocorrelation analysis (AC A) [26, 3lJ[32)- A 
full AC A over the BZ did not lead to a meaningful result, and 
hence we resorted to a restricted area ACA (red dashed lines 
in Fig. [3jd)) so as to extract the q x component. The results 
are shown in Fig. [3jf) and (g), indicating that q x = 0.0 has the 
maximum spectral weight. Thus, the experimental determina- 
tion of FSs gives us a SS vector of q = (27r/a)(0,0,0.86±0.1). 
While the MOKE results only provide the q z component, 
since q x = 0.0 from the present results, it indicates consis- 
tency between the MOKE and FS measurements within error 
bars for the obtained SS vector (Fig.[5Jh)). 

For the epitaxial ultrathin Fe/Cu(001) films, all available 
band structure calculations are thus inconsistent with the ex- 
perimental results, as they do not give nested regions. Also the 
magnitude of the SS does not match experiments. While theo- 
retical studies show that the SS in Erbium is due to FS nesting, 
the SS in 7-Fe is not expected to be derived from a single FS 
nesting vector [33]. The above results indicate the require- 
ment of an anisotropy between the k x -k y and k x -k z , which 
probably originates in the compressive strain along the c-axis 
l26l . Additionally, given that electron-electron correlations 
are not negligible in Fe (34), band structure calculations going 
beyond LSDA, e.g. LSDA + DMFT with SS order, are consid- 
ered necessary to explain the observations. Finally, a recent 
theoretical study has shown that the SS in fee Fe is a candi- 
date for spin transfer torque (STT: a current induced transfer 
of vector spin between magnetic layers) [33]. STT was pre- 
dicted for the incommensurate qi = (2tt / a)(0,0,0.6), but it was 
shown to be suppressed for q 2 = (27r/a)(0.5,0,1.0) due to the 
antiferromagnetic coupling between layers. Our results indi- 
cate an incommensurate SS vector of q = (27r/a)(0,0,~0.86), 
suggesting the epitaxial Fe/Cu(001) thin film is a candidate 
for STT. Since STT gives rise to fascinating properties: mi- 
crowave oscillations in a nanomagnet, current induced spin 
waves (35l [36), etc. and is technologically important for 
high- speed/high-density memory and switching applications, 
it would be interesting to investigate the same in ultrathin epi- 
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FIG. 3. In-plane and out-of-plane FSs of Fe(8 ML)/Cu(001) and 
autocorrelation analysis (a) and (b) In-plane and out-of-plane FSs 
measured over 1/4 of the BZ. (c) and (d) In-plane and out-of-plane 
FS crossings extracted from MDCs. FSs for full BZ were obtained 
by symmetrizing data from 1/4 of BZ. The regions indicated by the 
red dashed lines are the nesting parts and were used for the ACA. (e) 
MDCs for T-X along k x , k y and k z . The /c-cut positions are marked 
by black, red and green lines in (c) and (d). a=3.6\ and 3.55 A -1 are 
used for the in-plane and out-of-plane lattice parameters. The peaks 
observed only in the MDC along k z are indicated by the arrows and 
the corresponding EDCs are emphasized as thick lines in Fig. [TJf) 
and (g). (f) AC map for the restricted regions in out-of-plane FS. 
(g) Line profile of AC map at q z = (27r/a)x±0.86, indicated by red 
dashed lines in (f). (h) The comparison of the present result and early 
MOKE measurements: SS period 0ss = 2.3d (present result), 0ss = 
2.6d [16] and S s = 2.7 d fTSI , 
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taxial Fe/Cu(001) films. Furthermore, in-plane and out-of- 
plane FS mapping is expected to be important for studying 
strain control of properties in ultrathin films of elements, ox- 
ides, and their multilayers, hetero structures, interfaces, etc. 
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